Comparison of a new A resonance production model with electron and neutrino data 
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We consider resonance production by neutrinos focusing on the dominant resonance P33 at low 
energies with a detailed discussion of its form factors. The results are presented for free nucleon tar- 
gets. The A resonance is described by two form factors and and its differential cross sections 
are compared with experimental data. Further, we apply this approach to the electroproduction case 
and calculate its differential cross sections which are compared with electroproduction experimental 
data. Our approach to the analysis of A resonance is particularly simple and self-contained such 
that it could be helpful for the physical interpretation. 



I. INTRODUCTION 

The excitation of the resonances by electrons and neu- 
trinos has been studied extensively for a long time. The 
earlier articles Q, |M H, |^ tried to determine the 
pA transition form factors in terms of basic principles, 
like conserved vector current (CVC), partially conserved 
axial- vector current (PC AC), dispersion relations, etc. 
These and subsequent papers introduced dipole form fac- 
tors and in various cases other functional forms with 
additional kinematic factors in order to reproduce the 
data. As a result, the cross sections (differential and in- 
tegrated) were presented in terms of several parameters 
|6i5_7,_8,_9] . The relatively large number of parameters and 
the limited statistics of the experiments provided quali- 
tative comparisons but an accurate determination of the 
terms is still missing. A new generation of experiments 
is now under constructions aiming to measure proper- 
ties of the neutrino oscillations and they will provide the 
opportunity precise tests of the standard model. 

For these reasons it is important to improve the calcu- 
lation of the excitation of resonances with isospin 7 = 3/2 
and 7=1/2 looking into various terms that enter the cal- 
culations and trying to determine them, as accurately as 
possible. This has been done recently by three of us |10l |. 
In this contribution we will focus on the excitation of the 
A resonance which is dominant at the low energies con- 
sidered here ^3 ^'^'^ which is relevant for the study of 
single pion production in current and future long base- 
line (LBL) experiments, Hke K2k, MiniBoone, MINOS, 
CERN-GS, JHF, etc. We give expHcit formulas and dis- 
cuss form factors for P33(1232) where the formulas and 
form factors for other resonances Pii(1440), S'ii(1535) 
and 7)13(1520) can be found in 10]. 

In order to enlarge the event rates neutrino experi- 
ments use medium-heavy and heavy nuclei which brings 
in additional corrections such as the Pauli exclusion prin- 
ciple, Fermi motion and absorption and charge exchange 
of the produced pions in nuclei. However, in this contri- 
bution we restrict ourselves to free nucleon targets and 
refer to refs. for a discussion of nuclear effects in- 
cluded in our approach. 



The paper is organized as follows. In Sec. ^1 we give 
the general formalism for the production of the A res- 
onance following |lOj | emphasizing the minimal input, 
which is necessary. In Sec. IIIII we present numerical 
results of our formalism for neutrino production of the 
A++. In addition to the results presented in ^3 we 
compare with experimental results for the ratio of single 
pion resonance production (RES) and quasi elastic scat- 
tering (QE) which is important for understanding the 
problematic low region. Furthermore, we include sev- 
eral additional consistency checks of the theory by com- 
paring with more data for the spectrum of single pion 
neutrino production. In Sec. II VI we turn to a first com- 
parison of our simple approach with (fully differential) 
electroproduction data from JLAB and finally draw our 
conclusions in Sec. 



II. NEUTRINO PRODUCTION OF A 
RESONANCE: FORMALISM 

The double differential cross section for A resonance 
production is given by 



dcr 



dQ2dM^2 i67rA72 



i=i 



with Gp the Fermi constant and M the nucleon mass. 
The kinematic factors Ki{Q^ , E^,W) and the structure 
functions Wi (Q^ , W) which are expressed in terms of he- 
licity amplitudes are given in Ref. 0. The helicity am- 
plitudes depend on the Breit-Wigner factor 



fiW) 



(2) 



(W - Ma) - ^iTAiW) 

with T(W) = Fo qn{W)/q^{WR) and Fq = 120 GeV and 
form factors CY , ,i — 1 , . . . , 5 which are discussed in 
the following. 

For the matrix element of the vector current we use 
the general form 

75 u{p) f{W) (3) 
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C C C 1/ 

with Ax = -j^|-7A + -j^p'x + jj^Px-, C/ , i = 3, . . . , 6 the 
vector form factors, '0^ the Rarita-Schwinger wave func- 
tion of the A-resonance, and f{W) the s-wave Breit- 
Wigner resonance, given exphcitly in eq. iPJ. The con- 
servation of the vector current (CVC) gives {q^) = 
and the other form factors are determined from electro- 
production experiments, where the magnetic form factor 
dominates. This dominance leads to the conditions 



and 



Cr{q') = 0. (4) 



With these conditions the electroproduction data depend 
only on one vector form factor C^iq'^)- Precise electro- 
production data determined the form factor, which can 
be parameterized in various forms. 

Early articles describe the static theory and the 
quark model fis'l predicting the form factor for the ^NA 
vertex to be proportional to the isovector part of the nu- 
cleon form factors. Subsequent data [lj| showed that 
the form factor for A-electroproduction falls faster with 
increasing than the nucleon form factor which mo- 
tivated some authors to introduce other parameteriza- 
tions including exponentials and modified dipoles 
[l8l. The fmictional form 



Ml 



(1 



(5) 



4MS 



gives an accurate representation. In our approach [T3| 
we adopt this vector form factor and use CVC to deter- 
mine its contribution to the neutrino induced reactions. 
We can see in fig.^that our modified form factor in eq. 
(O agrees well with resent experimental results for the 
magnetic — A transition form factor G\j from Mainz, 
Bates, Bonn and JLab Details of the vector and 

axial contributions are presented in section lTTll where we 
shall estimate the contribution of from the electro- 
production data. 

The matrix element of the axial current has a similar 
parameterization 
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FIG. 1: The magnetic — A transition form factor GJ/ in de- 
pendence of Q2 16]. The solid and dotted lines denote our form 
factor eq. Q and the Sato-Lee(SL) form factor 17], respectively. 



modified dipole form 
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(1 



(7) 



i^f^Bxq^ - i^xq^Bf, 



+ ^^C^ + ^xq\^.C^ u{p)f{W) (6) 



with B\ = 



U^Px- 



The PCAC condition gives the relation C^(g^) = 
~M^1^ which for small = leads to the numerical 
value C^(0) = 1.2 0. The contribution of the form fac- 
tor Cq to the cross section is proportional to the lepton 
mass and will be ignored. 

The Q^-dependence of the form factors varies among 
the publications giving different cross sections and dif- 
ferent distributions even when the same axial vector 
mass Ma is used. For this dependence we shall use a 



The proton has a charge distribution refiected in the form 
factor. To build the resonance we must add a pion to the 
proton which creates a bound state with larger physi- 
cal extent. If the overlap of the wave functions has a 
larger mean-square-radius then the form factor will have 
a steeper dependence as is indicated by the electro- 
magnetic form factor for the excitation of A |Q| . Since 
the effect is geometrical we expect a similar behavior for 
the vector and axial vector form factors. For this reason 
we replace another factor used in previous article by 
the modified dipole in eq. © with 3Af^ ^ AMy. For 
the other two form factors C^{Q^) and Cf{Q^) we shall 



use = and CfiQ') 



\ct i. 



It is evident that 



there is still arbitrariness in the form factors with 
and being small. 
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III. NUMERICAL RESULTS 

We show in figure|21the relative importance of the var- 
ious form factors, where and dominate the cross 
section. The cross section from the axial form factors 
has a peak at = 0, while the cross section from 
turns to zero. The zero from the vector form factor is un- 
derstood, because in the configuration where the muon 
is parallel to the neutrino, the leptonic current is pro- 
portional to and takes the divergence of the vector 
current, which vanishes by CVC. The contributions from 
CY and are very small as shown in the figure |21 For 
this reason the excitation of the A resonance, to the ac- 
curacy of present experiments, is well described by two 
form factors and . 
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FIG. 2: Individual contributions to the cross section separated by 
form factors. The solid and dotted curves are the contributions 
proportional to the vector form factors C3 and C4 , respectively. 
The dashed and dot-dashed curves denote the parts proportional 
to the axial vector form factors and C4. 

An estimate of the vector contribution is also possible 
using electroproduction data. There are precise data for 
the electroproduction of the A and other resonances 0| 
including their decays to various pion-nucleon modes. In 
the data of Galster et al. cross sections for the channels 
(p -|- TT*') and (n -I- tt"*") are tabulated from which we con- 
clude that both / = 3/2 and 7=1/2 amplitudes are 
present. For instance, for W = 1.232 GeV the 7=1/2 
background is 10% of the cross section. 

For our comparison in ref. |lClj| we took the electro- 
production data after subtraction of the background, as 
shown in fig. O and then used CVC to obtain the con- 
tribution of V,t to neutrino induced reactions. We used 
the formula 



dV 



3 Q4 



cm, I—l 



dQ^dW 



dQ^dW 



(8) 



to convert the observed [l 
the reactions e+p ^ e + 



cross sections for the sum of 


to the vector contribu- 



pn 



n TT 



+ 



tion in the reaction // + p ^ /i^ + p + 7r+ denoted in eq. 
© by y . The factor 3/8 originates from the Clebsch- 
Gordan coefficients relating the matrix elements of the 
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FIG. 3: Cross section dV^ /dQ^dW for electroproduction in the 
A resonance in comparison with Galster et al. data (J^. The 
solid, dotted and dot-dashed lines are obtained with CY = 2.05, 
CY ~ 2.0 and C3 = 2.1, respectively. The dashed line denotes 
the backgroud contribution. 



two channels in the electromagnetic case to the matrix 
element of the weak charged current. We used the data of 
Galster et al. [li| at = 0.35 GeV^ and subtracted the 
background as suggested by them. Then we converted 
the points to the vector contribution for the neutrino re- 
action according to eq. ijHll. In the same figure we show 
the neutrino cross section with C^(0) = 2.05 (solid), 
CY(fi) = 2.0 (dotted), C|'(0) = 2.1 (dot-dashed) and 
contribution of background (dashed) and all other form 
factors equal to zero. Before leaving this topic we men- 
tion that the analysis of the electroproduction data [T^ 
included a contribution from the £'i3(1520) resonance 
which was found to be small. 

For the axial form factor we use the form given in eq. 
((7|). However, it is advisable to keep an open mind to no- 
tice whether a modification will become necessary. With 
the method described here we have all parameters for the 
A-resonance. We may still change the couplings by a few 
percent and vary My and Ma- 
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FIG. 4: Q -spectrum of the process vp n^p-K in comparison 
with BNL data 20, 21]. The solid curve is with Pauli factor and 
dotted curve is free case. 
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There is still the distribution 01 to be accounted 
for. The data is from the Brookhaven experiment 0,1^ 
where the experimental group presented a histogram av- 
eraged over the neutrino flux and with an unspecified 
normalization. We used the form factors for the A res- 
onance as introduced above and included the correction 
from a Pauli factor in the simple Fermi gas model with 
Fermi momentum pp = 0.16 GeV p.0]. The result is the 
solid curve shown in figure 0] which is satisfactory. For 
the relative normalization, we normalized the area under 
the theoretical curve for Q2 > 0.2 GeV^ to the corre- 
sponding area under the data. For the other parameters 
we choose 



cr(o) 



^3 

Mv 



1.95 
0.84 GeV 



C5^(0)-1.2, 

Ma = 1.05 GeV. (9) 
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FIG. 6: Ratio of distributions for RES and QE in comparison 
with BNL data 20]. The lines are the same as in fig.|5| 



In fact we made a x^-fit and obtained these values with a 
per degree of freedom of 1.76 for the complete re- 
gion. Furthermore, in order to reduce nuclear effects we 
performed a fit to all data with > 0.2 GeV'^ and the fit 
result gave a x^/d.o.f = 1.04. In the theoretical curves 
we av erag ed over the neutrino flux for the BNL experi- 
ment n%. The dotted curve is the calculation without 
Pauli factor and the solid one with Pauli suppression fac- 
tor included which has a small effect. It will be interest- 
ing to repeat this analysis when new data become avail- 
able. From fig. ^ we can see that in the region of small 
< 0.2 GeV^ the theoretical values are significantly 
above the experimental results which is not cured by 
a simple Pauli suppression factor which generates small 
corrections. Clearly, this discrepancy between data and 
theory is the origin of the worse jd.o.f. = 1.76 com- 
pared to the result of /d.o.f. = 1.04 where the low 
region has been excluded. 
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FIG. 5: Ratio of total cross sections for RES and QE in compari- 
son with BNL data 20] . The solid, dotted, and dashed lines have 
been obtained with = 1.08, = 1.04, and = 1.12, 
respectively. 

For the resolution of this problem it is reasonable to 
take the ratio of single pion production (RES) and quasi 
elastic scattering (QE) events, In figs. |S1 and [HI 



we show this ratio for total cross sections and distri- 
butions, respectively. The curves have been calculated 
with Ma — 1.05 GeV and various values for C^(0) and 
have been compared with BNL data [2(|. Since the ra- 
tio will reduce flux and experimental uncertainties it is 
an especially good test for a theoretical model. As can 
be seen from these figures we find very good agreement 
between our theoretical curves and the data, particularly 
also at low < 0.2 GeV^, because uncertainties in the 
low region which are common to both, resonance and 
quasi-elastic scattering, drop out in the ratios. 
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FIG. 7: The same as in Fig. fusing for the theoretical curves 
the PSY FFs 10] (red, solid and dashed lines) and SvH FFs J| 
(blue, dotted and dot-dashed lines) for various values of the axial 
mass Ma- 

We now turn to a more detailed comparison of our 
form factors (PSY FFs) jlfl] with form factors proposed 
by Schreiner and von Hippel (SvH FFs) and additional 
experimental data for the distribution of single pion 
production. Figs. |H1 and El show the flux averaged 
spectrum for A resonance production in comparison 
with BNL, FNAL, and CERN data [g^, IH IH Hi ■ The 
theoretical curves have been plotted using the PSY form 
factors (red, solid and dashed lines) and SvH form factors 
(blue, dotted and dot-dashed lines) for two Ma values. 
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FIG. 8: Same as in figure|7| with FNAL data 23]. 
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FIG. 9: Same as in figure|7lwith CERN data 



All these experimental results have been analyzed with 
a cut on the invariant mass, W < lA GeV, such that 
the A resonance dominates. One can see that the PSY 
FFs lead to a better description of the data than the SvH 
FFs for > 1 GeV^. Also notice that the form factors 
become insensitive to Ma for large where (which 
does not depend on Ma) is dominant. 



IV. A RESONANCE ELECTROPRODUCTION 

In this section we briefly discuss A resonance electro- 
production. We start with the fully differential cross sec- 
tion for neutrino reactions 5] as follows 



density matrix elements and Y™{9Tr, (p-n) spherical har- 
monic functions. This cross section can be easily con- 
verted to the electroproduction case by dropping axial 
vector parts and an appropriate change of the normal- 
ization factor. After some algebra we arrive at the fol- 
lowing expression for the electroproduction which has a 
part proportional to cos 2(f)* and cos 4>* and a part inde- 
pendent of (j)* 

- 2\/3 sin 9* cos 9* cos <j)* {K4D. 

- V^sin"^ 9*cos2(j)*{KeDe) 



2 

(11) 



with the normalization factor N 



9t. and 



the pion polar and azimuthal angles, respectively. The 
kinematic factors Ki, i = 1, 6 and the structure func- 
tions Wi , Di can be found in and depend on the vector 
form factors introduced in Sec.lTTI 
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FIG. 10: Electroproduction data 25] for the reaction e + p 

e+p + n". 



dV 



1 



dQMW^dn* 



4^dQ2dM^2 



J- n 



2 _ 1 

P(3 _i)Rey2' 



P(3i)Rey2^ 
(10) 



,4 

with the double differential cross section ^g" ~ 
^Y.^t=iKiWi given in eq. (|TJ, ^(33), p(3i), p(3 _i) spin 



For our numerical comparison with electroproduction 
data we have calculated the fully differential cross section 
in eq. using the vector form factors given in eqs. 
0) and ||SJ). In fig. ^| we show electroproduction data 
125] for the reaction e + p ^ e + p ~\- n'^ in dependence 
of the pion polar angle 9^^ for various pion azimuthal 
angles (j)^ for fixed = 2.8 GeV^ W = 1.235 GeV 
and electron energy — 3.2 GeV. Similarly, in fig. ITTI 
experimental results at small momentum transfers = 
0.45 GeV2 and Q'^ = 0.75 GeV2 are shown. As can be 
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FIG. 11: The same as in fig. [THIfor small data 126 



seen the description of the data by our simple model is 
already quite reasonable with exception of data points 
at |07r| > 135° in the upper panels in fig. E| It should 
be noted that the theoretical curves have been obtained 
without fine tuning of the parameter 6*3^ (0), for which we 
took the value C^(0) — 2.0, and without including the 
background from the 1=1/2 channel and a non-resonant 
background which is expected to be of the order of 10%. 
A more detailed comparison will be presented elsewhere 



There exist already more complicated calculations, i.e., 
with more parameters and a more complicated functional 
form of the form factors. On the other hand, our ap- 
proach is particularly simple and self-contained. There- 
fore it is helpful for the physical interpretation. 

We have analyzed the differential cross section in 
terms of two form factors which are described by two 
free parameters Ma and My and compared our results 
with several experimental data. We obtained a good fit to 
BNL data ^ with a x^/d-o.f = 1.04 for > o.2 GeV^ 
However, for < 0.2 GeV'^ the theory value is too high. 
A possible explanation are medium effects which are not 
properly taken into account. The Pauli suppression fac- 
tor has a small effect which does not resolve the discrep- 
ancy. To shed light on this problem we took the ratio 
'ct(X^? ^^'^ found very good agreement with BNL data 

m 

From these comparisons we conclude that the neutrino 
production of the A resonance can be described by two 



form factors 



and C^. 



production data and C5 
neutrino data. 



is consistent with electro- 
can be determined only from 



We verified the consistency of our model by comparing 
with further existing neutrino production data. More- 
over, we found that our form factors gave a better de- 
scription than the form factors originally proposed by 
Schreiner and von Hippel [|. 

Finally, we performed a first comparison with more 
differential electroproduction data and found reasonable 
agreement also here. In this comparison we have ne- 
glected the background from the isospin 1=1/2 channel 
and a non resonant background which is of the order of 
10% atW = 1232 GeV. 



V. CONCLUSIONS 

Single pion neutrino production at low energies is im- 
portant to measure properties of neutrino oscillations. 



In the future we will include higher resonances and 
non-resonant background and perform a more detailed 
comparison with electroproduction data in order to study 
further the connection between neutrino- and electropro- 
duction. 
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